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Le domalne d'application de I'invention est celui de la prospection sismique. 
L'invention conceme rimagerie sismique, et plus partlculierement I'estimation du 
taux de couverture d'illumination dans le domalne 3D migr6 necessaire pour 
I'obtention d'une image haute resolution, ^ amplitude r6eHe, de la structure du 
5 sous-sol. 

La prospection sismique consiste d'une maniere g6n6rale, ^ 6mettre dans le 
sous-sol, a I'aide d'une ou plusieurs sources sismiques, des ondes sismiques, d 
enregistrer en surface, en fonction du temps, des donn6es sismiques 
correspondant aux ondes sismiques r6flechies sur les Interfaces g6ologiques du 
10 sous-sol (encore appel6es rSflecteurs) ^ I'aide de r^cepteurs (encore appetes 
g§ophones ou hydrophones selon que Ton prospecte ^ terre ou en mer) puis a 
tralter ces donnees pour en extraire des informations utiles quant ^ la g§ologie du 
sous-^ol. 

Le traitement des donn§es (encore appel6es traces sismiques) enregistr^es 
15 par les r^cepteurs pemriet en particulier de concevoir des images presentant des 
informations relatives d la g^ologie du sous-sol. / 

L'imagerie sismique a haute resolution, a amplitude r^elle, a pour objectif de 
fournir une image du sous-sol qui est non perturbee notamment par les effets de 
la propagation du champ d'onde dans le sous-sol et par les effets de I'acquisition 
20 (on parle alors d'empreintes d'acquisition). 

La representation en image de la r6flectlvit§ acoustlque du sous-sol est 
r6all8§e par un proc§d§ dit de migration qui pemriet de restituer les forrtie^ 
correctes des interfaces gdologiques. Les donn6es sismiques sont pour cela 
converties depuis le domalne dans lequel elles sont acqulses vers le domalne 
25 migr6 dans lequel les Interfaces g6ologiques sont representees S leurs veritables 
positions. 

Dans le domalne mlgr6, rechelle verticale peut etre le temps (on parie alors 
de migration temps ou « time migration » selon la terminologie anglo-saxonne) ou 
la profondeur (on parle alors de migration profondeurou « depth migration » selon 
30 la terminologie anglo-saxonne). 
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_ ®^ gen^ralement r6allsee avant sommatlon (on parle ajprs cie. 

migration avant somme ou « pre-stack migration » seion la termlnologie angle-' 
saxonne). avec un mappage direct du signal sismlque depuis les traces slsrrilques 
vers le domalne migre. 

Udentiflcation precise de la g^ologie du sous^ en un point particulter est 
dimotement re,«e 4 la capacity de remission sismique a illuminer ce point Si iin 
po,nt du sous-sol est hautement iilumln*. le «lev6 sismique permet d'^valuer en 
toute vraisemblance la g^,o0e du sous-sol en ce point. En .evanche, si un point 
du sous-sol est peu Ulumln.. 11 est possible que le relev. sismique fbumlsse une 
Evaluation imprMse. voire Inconecte, de la gSologle du sous^ en ce point 

Afln d'obtenlr une .eprSsentalion p,«cise du sous-sol. les amplitudes r^elles 
S..^!!tt" '"^^ parfaltemen. 

d-amlt """^ ^"^^""^^ ^' *angement 

damplrtude en un potnt particulier du sous-sol est le rfeultat d'un changement 

parlicuiiere en ce point 

Mais I'lllumination des nMlecteurs est g^n^r^iement non unilbm,e. Ble 

depend en particulier : 

- de la configuration de Tacqulsltion des donn6es sismlques • 

- des variations laterales et vertlcales de la vltesse de propagation des ondes 
sismiques dans le sous-sol. 

Afin d'obtenir une in^ l,aute r&olution de la r^iecBvltE du sous-sol 11 est 

-^"--^ ''•effecue; une 

™on d amplitude pour prendre en compte la non-unfcrmit. d il,uminat,on des 

I'aide^Tr" ' Pr^cedemment. cette ocnection (,«alis^ ^ 

I aide d un fecteur de pondSratlon) est n^cessaire car • 

- pour Imager pr^cis^ment la reflectivity du sous-sol, 11 faut corriger les 
d,stors,ons en amplitude li.es , la propagation du champ d'onde dans le sous- 
sol (divergence sph6rique, absorption, directivity, etc.) : 
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- pour am§llorer la resolution de I'lmage de r6flectivit6. il faut prendre en compte 
la variabilite du taux de couverture d'lllumination de migration pour chacune des 
IncHnaisons d'un r§flecteur qui sont prises en compte lors de la migration. 

Le taux de couverture d'illuminatlon correspond, pour un point image donn6 
5 du domaine migr§, au nombre de paires (source, r^cepteur) dont les contributions 
pour la constitution de I'amplitude de I'image de r^flectivite Interfdrent 
constmctivement (condition de phase stetionnaire). 
II a §t6 propose dans le document : 
« Bleistein, N. 1987, On the Imaging of Reflectors In the Earth: Geophysics, 
10 Soc. Of Expl. Geophysics, 931-942 », 

une methode permettant d'obtenir, pour une o^nfiguration reguli§re 
d'acqulsition, un facteur de pondSration d'amplitude permettant de rSaliser une 
representation haute resolution, en amplitude r6elle de la reflectivite du sous- 
sol. 

15 Cette methode est une methode « inverse » constituee d'une etape de 

modeiisatlon (modeiisation de traces sismiques d partir d'un module de 
reflectivite), suivie d'une etape de migration (depuis les traces sismiques 
modelisees, vers une reflectivite imagee), le facteur de ponderation etant alprs 
choisi de maniere a obtenir une concordanqe parfaite entre la reflectivite imagee 

20 et le modeie de reflectivite. 

En suivant cette « methode inverse » proposee par Bleistein, et comme cela 
est represente sur la figure 1 , 11 apparatt que : 

- pour une configuration d'acqulsition donnee (repartition des paires (source S, 
recepteur R)), 

25 - pour un point Image x donne dans le domaine migre, 

- pour une inclinaison de reflecteur donn6e (indiquee par un vecteur p de 
pendage), 

le taux de couverture d'illuminatlon l(jc,p) est determine par le nombre de paires 
(source, recepteur) tel que leurs contributions de migration interferent 
30 constmctivement ^ I'image du reflecteur de pendage p au point image x . 



ler cl6p6t 



4 



Les expressions d6veloppees par Blelstein pour le taux de couverture 

d'lllumlnation \Cx,p) integrant en une seule fois reffet de la totallte de la 
configuration g§om§trique de I'acquisition. 

Par consequent, la m6thode propos6e par Bleistein ne peut §tre utilisee que 
pour une distribution regull^re et param6trique de I'ensembie des palres (source 
recepteur). 

Cette condition restrictive d'application en limfte par consequent futilisation a 
un nombre redult de configurations d'acquisition. 

Pour des configurations d'acquisition a geometries irregulieres. II est done 
necessaire de definirles taux de couverture d'lllumlnatlon individuels l(i JJ. r), 
pour chacune des palres (source s, recepteur ?) impliquees dans la migration' 
Lesdits taux individuels contribuent bien entendu au taux de couverture 
d'illumination total I(^ ,;)selon I(x,p) = j;^i(x,p;s,P)). 

L'^de condulte par Bleistein est toutefbfe 4 la l,ase cTun nombre de 
developpements qui ont eu pour but d'estimer le taux de couverture d'lllumlnalion 
dans te domalne mIgrS. pour n'importe quelle configuration g4om6trique de 
I acquisition, et pour nlmporte quelle distribution des paires (source. ,*cepteur) 

Au cours des ann^es. il a ^ propose diverses appmximaUons pour l onde 
sismique. qui est raiment ot« pour Testimation de la condition de phase 
stationnaire (condition d'interfirence constructive). .^ 

Les premieis dSveloi^ments ont consists 4 modaiser par lancer de rayons 
(« ray tracing » salon la tem,inologie anglo-saxonne). pour ct,acune des paires 
(source r*:epfeur 'r). les points de rSflexlon spficuiaires te long d'un horizon 
(carte 2D) sSlectionn^ dans te domains migi6 (3D). 

Cette 6tape de modaisation de la r^flexton permet de d^temiiner les 
quantity sulvanles : 

- La localisation des points de rMexIon x.CsTr) te long de fhortzon 
selectionne ; 

- Le temps de trajet tr(ir ■.s,r)6e reflexion speculaire ; 
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Le temps td de trajet de diffraction le long de i'liorizon selectionn6 peut alors 
§tre estim§ par d§veloppement en s6rie de Taylor du temps de trajet tr de reflexion 
autour du point xr(s,r) de reflexion sp6culaire (d6tennination des d6rlvees 
premiere et seconde VjX^r) ®* ^xjXxr))- 
5 La figure 2 lllustre cette procedure au cx)urs de laquelle le temps de trajet de 

diffraction td(x) au point image x est estime a I'aide du temps de trajet de 
reflexion tr(xr) et plus particulidrement par son d§veloppement en serie de Taylor 

autour du point de reflexion sp6culaire xr . 

La sommation de tous les mappages d'ondelette obtenus pour chacune des 
10 palres (source, r§cepteur) pemiet ensuite d'estimer le taux de couverture 
d'illumination dans le domaine migr§, le long de I'horizon selectionn6. 

En mod6lisant I'ondelette sismique selon un pic impulsionnel, d bande de 
frequence infinie, il a et6 montr§ dans le document : 

« Cain. G., Cambois, 6., Gehin, M. et Hall, R., 1998. Reducing risk in seispiic 
is acquisition and Interpretation of complex targets using a Gocad-based 3-D 

modeling tool, 68th Ann. Intemat. Mtg: Soc. of Expl. Geophys., 2072-2075 ». 
qu'il est possible d'obtenir une carte d'Impact, seule la contribution speculaire ^ur 
un liorizon s^lectionne 6tant alors effectivement prise en compte. 

En prenant en compte une approximation realiste d'une ondelette sismique a 
20 bande limitee, il a §te montr6 dans les documents : 

- « Schneider, W., Jr^nd Winbow, G., 1999, Efficient and accurate modeling 
of 3D seismic illumination, 69th Ann. Internat. Mtg: Soc. Of Expl. Geophys, 
Expanded Abstract. 633-636 » ; 

- « Laurain, R. et Vinje, V., 2001, Pre-stack migration and illumination maps, 
25 71th Ann. Internat. Mtg: Soc. Of Expl. Geophys, Expanded Abstract, 929- 

932 » ; 

qu'il est possible d'am6liorer I'estlmatlon du taux de couverture d'illumination le 
long d'horlzons s§lectlonn6s. 
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particulier sur un balayage d'horizons 

s6lectlonn6s. se sont montrees relatlvement fiables et precises. Leur efficacite est 
cependant Iimlt6e, notamment du fait que : 

• la compensation du taux de couverture d'llluminatlon est Ilmlt6e & un nombre 
fini d'horizons selectionnes et I'approche. au contraire de celle d§veloppee par. 
Bleistein, n'est done pas totalement 3D ; 

• pour chacune des palres (source I, r^cepteur ?) impliqu6es dans ie 
migration, une 6tape de modellsation d§dl§e dolt §tre reafis^e apr^s la 
migration. Or cette 6tape de mod6llsatlon, par lancer de rayons, pour la 
determination des quantit^s (3c,,V,.,,A,y,) mentionn6es pr^cedemment. Ie 

long de chacun des horizons selectionnes. est partlculidrement coQteuse en 
temps. 

Afin de surmonter ces limitations, et de se rapprocher d'une apprache 
volumique 3D. une solution, en deux etapes. pour I'estlmatlon puis la 
compensation du taux de couverture d'llluminatlon a 6t§ propos^e. en particulier 
dans les documents : 

- « Albertln. U.. Bloor. R.. Beasley. C. Chang. W.. Jaramilfo, H.. Mobley. E. 
and YIngst. D., 1999. Aspects of true amplitude migration. 69th Ann. Intemat. 
Mtg: Soc. Of Expl. Geophys. Expanded Abstract. 1247-1260 » ; 

- « Audebert. F.. Froldevaux, P. Huard. I.. Nicoletis. L. and slay-Lucas. J., 
2000. A mutll angle toolbox for restored amplitude images and AVA gathW 
70th Ann. Intemat. Mtg: Soc. Of Expl. Geophys. Expanded Abstract 1020-' 
1023 » . 

La premiere §tape consiste. dans Ie domaine mlgr6. en une analyse du taux 
de couverture d'illumlnatlon multi-pendage utillsant les tables pr6-etablles pour la 
migration de KIrchhoff des temps de trajet de diffraction Uix-.s ,r). 

Lors de la seconde etape. les Inverses des taux de couverture d'illumlnatlon 
estimes sont utilises afin de realiser la sommation ponderee des donn^es 
sismiques. 

Pour cheque point Image x et pour toute une gamme de pendage p de 
reflecteur. il a ete propose d'Incrementer un compteur d'occurrence d'illumlnatlon 
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I(x^p) , chaque fois qu'une contribution de migration attaint le point image x avec 
ia composante p de pendage. 

De mani§re similaire S ce qui avait 6t6 propose par Cain et al (cf. ci-dessus), 
le critere utilise pour augmenter le compteur d'occurrence d'illumination prend en 
5 compte un modele impulsionnel, a bande de frequence infinie, de Tondelette 
sismique. 

Cette approximation hautes frequences est telle que, pour une composante 
de pendage donn§e, seules les contributions de migration speculaires 
correspondant a ladite composante de pendage sont prises en compte pour 
10 Testimation de Illumination. 

Et cette approximation hautes frequences pr6sente un certain nombre de 
limitations : 

- Testimation de rillumination est bials6e puisque la nature a bande limitee du 
signal sismique n'est pas prise en compte**; 
15 - Afin d*6viter la division par des zeros, les taux de couverture dlllumination 
sont ranges par classes de pendage. La taille des classes est toutefois un 
element critique. Si la taille de la classe est trop petite, les risques 
d'instabilite numerique sont plus grands ; si la taille de la classe est trop 
grande, la precision quant a Testimation de ia couverture dllluminatioR est 
20 diminuSe, > 

L'invention vise un precede d'estimation du taux de couverture d'illumination 
sismique dans le domaine 3D migr§ non sujet ^ ces limitations. 

Plus pr6cisement,ifrinvention a pour objectif de proposer un proc6d6 qui 
prend en compte la nature a bande Iimit6e du signal sismique afin d'estimer des 
25 taux de couverture d'illumination plus realistes, et qui. dans le meme temps, 
resout le probldme de la tailte des classes. 

A cet effet, Tinvention propose, selon un premier aspect, un procede 
d'estimatlon du taux \{x^p ) de couverture d'illumination sismique dans le domaine 

3D migre en un point image x, pour un pendage de vecteur ^ , caracteris6 en ce 

30 que, pour chaque paire (source , recepteur r ) du releve sismique, il 6value le 
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de couverture dNIIurnlnation rejatif ^ . ladlte paire. (source J, 
r6cepteur r ) par mise en oeuvre des 6tapes de : 

• d6teimlr.ation du temps de trajet de r6flexion t,( (J ) ; ^ . ; ) depuis la source 
3 , jusqu'au point de reflexion sp6culaire sur le reflecteur plan passant par ie 
point image x et perpendlculaire au vecteur ? de pendage. puis de retour vers 
le reflecteur r , 

a partir du temps de trajet de diffraction ^x; s, r) depuis la source 2 
jusqu'audit point Image x, puis de retour jusqu'au reflecteur 7 ; 
. d'Incrementation dudit taux I(x,p',s,r) de couverture d'illumlnation relatif a 
ladlte paIre (source 2 , r^cepteur ?) en fonction de la difference entre le temps 
de trajet de diffraction ta(i; 2. i^) et le temps de trajet de reflexion tr(x, (J ); J. 

Des aspects preferes. mals non llmitafrfs du precede selon le premier 
aspect de I'Invention sont les suivants : 

- le precede comporte en outre une etape de sommation de chacun des taux 
Hx,p;s,r) de couverture d'illumination relatif a une paire (source s, recepteur 

r). de maniere d determiner le taux de couverture d'illumlnation total 
Hx,p) = J^Iix,p;s,r). 

- au cours de retape< d'Incrementation. le taux Iix,p;s,f) de couverture 
d'illumlnation peut etre Incremente en utillsant une fonction d'Incrementation 
i(td. tr; s,r), selon mp;s,r) = I(x,p;s,P) + i(td, f,; ?,F), 

- la fonction d'Incrementation I peut etre fonction de I'ondelette sismique s(t), 

^ la fonction d'incrementatron I peut s'exprlmer en fonction de la derivee selon le 
temps de I'ondelette sismique, 

- II est envisage de tenir compte d'une correction a priori w(x;s,r) du taux de 

couverture d'illumination. correction qui est generalement mise en jeu loi^ de 
la migration, 



^ r6tape de d6termlnation comprend le d^veloppement en s6rie de Taylor, du 
second ordre, du temps de trajet de diffraction U{x; s, r) autour du point 
Image x , 

- le point de reflexion specuiaire x,(p) est detennin^ le long dudit r6flecteur de 
sorte que le temps de trajet de diffraction audit point reflexion specuiaire x^(p) 
est stationnaire, 

- I'etape de determination utilise les cartes des Isochrones de migration U{x', s, 
r ) specifi^es pour toute paire (source j , r6cepteur r ) impliqu§e dans la 

migration et tout point image jc du domaine 3D migr6. 

- I'estlmation du taux \{x,p) de couverture d'illumination sismique dans le 
domaine 3D migr§ est realist lore de la migration par sommation de Kirciilioff 
des donnees sismiques enregistrees au coure du releve sismique 3D. 

Selon un autre aspect, I'invention propose un proc6de de correction des 
amplitudes des donn6es sismiques enregistrees au cours d'une prospedion 
sismique 3D afin de compenser I'effet d'une illumination non uniforme des 
reflecteurs du sous-sol, caract6rise en ce qu'll comporte les etapes consistant a .: 

• estimer le taux de couverture d'illumination \{x,p) par le precede selon le 

premier aspect de I'invention, 
. utiliser I'inverse \'^{x,p ) dudit taux comme facteur de pond6ration a appliquer 

d chacune desdites amplitudes des donn§es sismiques. 

Selon encore un ^utre aspect, I'invention propose un proc§d§ pour la 
selection, en fonction de la cible d'une prospection sismique 3D, d'une g§om6trie 
d'acqulsition parmi une pluralite de g§om§tries d'acquisition, caract6ris6 en ce 
qu'il comporte les etapes consistant a : 

• determiner, pour chacune des geometries d'acquisition considerees, le taux 
de couverture d'illumination \{x,p ) par le precede selon le premier aspect de 
I'invention, 

• s6lectionner la g6ometrie d'acquisition procurant le taux de couverture 
d'illumination optimum en fonction de la cible. 
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- .^.^.P?^*^.'.^"*.s et avantages de I'invention apparaltront mieux a la 
lecture de ia description detalHee suivante d'une forme de realisation pr6fer6e de 
celle-ci. donn§e ^ titre d'exemple non llmltatif et faite en r6f6rence aux dessins 
annexes sur lesquels : 

- la figure 1. deja commentee, illustre la contribution constructive ou non, en 
fonction de s(td - 1,), de la trace sismlque migree S {'image du r6flecteur peht§ 
en X \ 

- la figure 2. d6jd commentee, illustre I'estimatlon r6alis6e dans le cadre de 
r6tat de la technique du temps de trajet de diffraction au point image U^x) a 
partir d'une developpement en serie de Taylor du temps de trajet de reflexion 
au point de reflexion sp6culaire tr(xr) ; 

- la figure 3 illustre le proc6d6 selon Tinvention selon lequel on estlme le point 
de reflexion sp6culaire Xr ainsi que le temps de trajet de reflexion t^^,) k 
partir d'un developpement en s6rie de Taylor du temps de trajet de diffraction 
au point image U(x)- 

- La figure 4 est un organigramme pr6sentant de manldre sch6matique les 
principales etapes du proc6d§ selon I'invention. 

L'invention propose un proc§d6 d'estimation du taux de couverture 
d'lllumination dans le domalne migr§ /(5,^;5,F) pour n'importe quelle position 
mlgr6e x, n'importe quelle composante de pondage p et n'importe quelle 
distribution des traces slgmiques. en utilisant uniquement les quantlt6s les plus 
basiques mises en jeu dans la migration par sommation de KirchhofF. k savoir : 
• les coordonn^es (^.P) des paires (source, recepteur). I representant la 

position d'une source S et r celle d'un r§cepteur R ; 
. les cartes des Isochrones de migration t. ; ^ . P ) pour chacune des paires 
(source s , recepteur r) impliquees dans la migration. 

L'invention s'appllque sans distinction k I'imagerie en temps et h I'lmagerle 
en profondeur. 
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L'invention propose un proc§d6 du type comprenant une etape de 
mod^lisation suivie d'une etape de migration. 

En ce qui conceme r§tape de mod§lisation du point de reflexion xr et du 
temps de trajet de reflexion i^xr), le proced6 objet de l'invention applique la 
condition de phase stationnaire a un developpement en serie de Taylor des temps 

de trajet de diffraction U{x). 

L'ordre dudit developpement est un developpement du second ordre. 

En ce qui concerne l'§tape de migration, le proc6de objet de l'invention utilise 
les cartes des isochrones de migration disponibles pour la migration 
classiquement r^alisee par sommatlon de Kirchhoff. 

Pour chacune des palres (source 5 . r6cepteur r ), le taux de couverture 
d'illumination l(3c,p ,1 ,r) est mis a Jour, en reference a la figure 3, selon I'^quation 
suivante : 

i(3c,^;^.r) = l(5,^;s.;) +i(td,tr; J,r) Equation(1) 

Avec 

• 3c , le point image dans le domaine migr§ x = (xi, X2,X3) ^ ; 

• p le vecteur de pendage repr6sentatif de la composante de pendage 

selectionn6e, p = (pi, P2,P3) ^ ; • 
, les coordonnees de la source, s = (si, sz, S3) ^ ; 

• r, les coordonnees d(rr6cepteur, r = (ri, ra, xz) ^ ; 

. U (x ; s , r), \e temps de trajet de diffraction depuis la source s, jusqu'a la 
position migr§e x, etde retour jusqu'au recepteur r (trajet indiqu§ par td sur la 
figure 3 et repr6sente en trait plein) ; 

• tr(xr(J ) ; s, r), le temps de trajet de reflexion depuis la source s, jusqu'au 
point Xr de reflexion sp^culaire sur le r6f!ecteur plan passant par le point 
Image x et perpendiculaire au vecteur J de pendage, et de retour jusqu'au 
recepteur r (^jet indiquS par tr sur la figure 3 et repr6sent6 en trait polntille) 
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• i(td, tr ; s,r ) est une fonction d'incrementation dpnt rexpresslon sera d6taill6e 
ci-apres. 

Comme cela a d6ja 6t6 mentlonn6 pr6cedemment. le proced6 selon 
nnventlon utilise avantageusement ies cartes des Isochrones de migration t^ {x ; 
s, r) mises en jeu dans la migration par sommation de Kirchhoff. pour chacune 
des paires (source s, r^cepteur O Impliqu6es dans la migration. Ces cartes sont 
ainsi des « donn§es d'entree » du precede selon I'invention, 

II utilise 6galement Ies d§riv6es premiere et seconde V,^,(x;?.r). 
A,j,(x;s,r) des isochrones de migration lorsque celles-ci sont disponibles. Dans 
le cas contraire. lesdites derivees ne sont pas des donnees d'entree. mais 
peuvent etre calculees dans le cadre du precede selon invention 

La figure 4 est un organigramme repr6sentant Ies diff§rentes §tapes du 
proc6d6 selon un mode de realisation pref§r6 de I'invention. 

Ledit proc6de consiste S r6aliser. pour chacune des palres (source 
r6cepteur impliqu^es dans la migration, une 6tape 1 mise en cBuvre pour toui 
point image x du domaine 3D migr§. 

Ladite etape 1 comprend, pour toute composante de pendage p 
s6lectionnee pour le calcul de I'illumination I(x,p;s,r) , deux operations 2. 3. 

La premiere operation 2 consiste S determiner le temps de trajet e^e 
reflexion \r{xrCp); s. ?).pu point de reflexion speculaire x,(J) sur le r6flecteur 
plan passant par le point image x et perpendiculaire au vecteur p de pendage. 

La determination de I'operatlon 2 est realisee a partir du temps de trajet de 
diffraction U{x;s,r) connu des cartes des Isochrones de migration. 

La seconde operation 3 consiste ^ Incrementer le taux de couverlure 
d'.lluminatlon dans le domaine migre I(x,p;s,r) d'une fonction dWmentaflon 
i(td. tr ; J,F ). ladite fonction d'incrementation prenant en compter 

- la difference { U (x ; s, 'r) - t^CMp) ; s, ?) }.entre Ies temps de trajet 
mentionnes cl-dessus ; 
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- Tondelette sismique s(t), 

^ et, 6ventuellement. le terme d^amplitude de la migration w(x;s,r) . 

A Tissue des operations 2 et 3, le taux d'amplitude d'illumination T(x,p;SJ) 

est connu pour un point image x donn6 et une composante p de pendage. 

En repetant les operations mentionnees ci-dessus pour chacune des patres 
(source s, r6cepteur r), on determine finalement le taux d'illumination dans le 
domaine migr6 \{x,p ). 

On additionne pour cela Tillumination relative a chacune des paires (source 
s , recepteur r ) selon = ^/(x,^;?,F) . 

Le point de reflexion Xr , ainsi que le temps de trajet de r6flexion tr sont 
d6termines (au cours de I'operation 2 mentionn6e ci-dessus) a partir d'une 
analyse en phase stationnaire, par un developpement en s6rie de Taylor, du 
second ordre, du temps de trajet de diffraction td autour du point image x . 

A partir de cette analyse (cf. Equation (2) ci-dessous), on obtient les 

expressions suivantes pour Xr ettr (equations (3) et (4)) : 

tAx;SJ) = Ux/J,fO+(^jJix-x,)+^ix-x,f(^ Equation (2) 



Xiip)-x -M, . b Equation (3) 



U^KiP)' ^. r) = i6{x : s, r)-i.6^F-^6 Equation (4) 

00 dans les Equations (3) et (4), 

• M est une matrice (3 x 2) d§crite par deux vecteurs s'6tendant le long du plan 
de reflexion. Ces deux vecteurs sont par consequent perpendiculaires au 
vecteur de pondage p. 

• b est un vecteur (2x1) des deriv6es du second ordre du temps de trajet de 
diffraction le long du plan de reflexion, 

6 = M ^ ) Equation (5) 
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. F est une matrlce (2 x 2) des derives du premier ordre du temps de trajet..de_. 
diffraction le long du plan de reflexion, 

F = M^(A,^/,).M Equation (6) 

On decrit cl-apr§s I'obtention des Equations (3) et (4). 
5 On precise id que le vecteur ~p de pondage represente le pendage apparent 

d'un r6flecteur dans le domaine migr6. ce r§flecteur passant par le point image x . 
Le vecteur p de pendage s'exprime ainsi selon : 

P = (ar3/ax, .3^3/3x3 ,1) Equation (7a) 

QKjdx, represente le pendage apparent le long de I'axe xi et Qx,ldx^ , le pendage 
10 apparent selon I'axe Xa. 

0 et V/ represented respectivement I'incllnalson et I'azimut du r^flecteur 
pent§ passant par le point image x ou I'iHumination est 6valuee. 
On a ainsi les expressions suivantes: 

COSe = (l + p^+P2)-^'2 

Equation (7b) 

15 tg¥=p2/Pi 

Equation (7c) 

A partir du d^veloppement en serie de Taylor du temps de trajet de 
diffraction (connu de la carte des isochrones de migration), on a. comme on I'a vu 
precedemment : 

1<l(3^.) = t.(3c) + (V,^,)-(^,-3c)■^i(^,.J^)^A,^.,.(■c,-^)-^... Equation (2) ^. 

Le point de reflexion- x, est tout d'abord estim§ le long du plan de reflexion 
pour lequel le temps de trajet de diffraction est stationnaire (condition pour une 
reflexion sp§culaire). 

Afin de ne consid6rer que les positions appartenant au plan pent6. on 
introduit un nouveau vecteur x , tel que : 



20 



25 x^=x+MJCp=x + 



^ 1 0 ^ 
0 1 

-Pt -P2J 



Equation (7d) 
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On note Ici que, comme d6sir6. les deux vecteurs de la matrlce M sent bien 
perpendlculaires au vecteur p de pendage. 

La re6criture de I'equation (2) en utilisant I'^quation (7d) fournit I'expression 
suivante : 

f.(3c,) = /,(J) + (v,«.)'M.x^ +l^J.M^A,^i,.M^;: +... Equation (7e) 

soft: f,(3c,) = «<,(3c)+5^Jc,+i4.F^;+... Equation(7f) 

Dans cette Equation (7f) : 

5^ ^^,t,jM ,d'oCi 5=M^(v,<J, et F^M\t^,j^M. 

On en d§duft par consequent la position 3c, pour laquelle le temps de trajet 
de diffraction f^jEp) ®st stationnaire : 

^Jt^i^A =0,=^^ =-F-'5 => x,=5 + MJc, Equation (7g) 

En combinant les equations (7g) et (2), on obtient I'expression du temps de 
trajet de reflexion de l'6quatlon (4) : t, =t^ix,')»t^{^x)-^b'^ JF'^b . 

Comme cela a Me expos6 precedemment, la taux de couvertcire 
d'illumination est Incr6mente, au cours de I'operation 3, en fonction de' la 
difference entre le temps td de trajet de diffraction et temps tr de trajet de reflexion. 

Meme si le temps de trajet de r6flexion tr(5^(p ) ; ? , r ) differe du temps d^ 
trajet de diffraction td(x ; "I , r). du fait de la nature a bande Iimit6e des donn^es 
sismiques, la migration de la reflexion ayant lieu a peut contribuer 

constructivement h I'image du r§flecteur pent6 en jc . 

La selection entre contributions de migration illumlnantes (interferences 
constructives) et contributions de rhigration non Illumlnantes (interferences 
destructives) est r§alisee grace la fonction d'incr§mentation i(td, tr; ?,F). 

Selon un mode de realisation de I'invention, la fonction d'incrementation 
s'exprime selon i(td,tr; J,F ) = i(/rf(x;J,r)-f^(x,(p);?,F)), oO : 
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- ^(0 est I'ondelette sismique, qui peut §tre. monpchrpmatlque. a bande.. 
limitee. ou a bande inflnie et qui peut etre cholsie d6pendante du temps et 
de I'espace. 

- m est la deriv6e selon Ie temps de I'ondelette sismique i(0 = d[sit)]/dt . 

5 Une d6riv6e temporelle est generalement appllquee a chacune des traces 

sismlques avant la migration de KIrchhoff (en temps ou en profondeur) AInsi 
notamment dans un souci d'unifomilt6 des quantlt6s etudi6es. on choisit la derivee 
m selon Ie temps de I'ondelette sismique. au lieu de I'utllisatlon de I'ondelette 
sismique ^(0 elle-meme. 

10 Selon un autre mode de realisation de I'invention. notamment loi^que les 

traces sismiques ne sont pas d6riv6es selon Ie temps, la fonctlon d'incr^mentatlon 
prend en compte I'ondelette sismique .(/) au lieu de la d§riv6e i(0 de celle<:i 
selonletemps;onaalorsi(t..tr; ?,F)= sitAx;s,n~UK(py,s,n). 

Un temie de correction d'amplitude de la migration MxJ,r) est 
15 classiquement calcule lors de la migration, et prend notamment en compte une 
g6om§trie reguliere d'acquisition. 

Ce temne y.(x;s,F) peut done §tre avantageusement utilise dans Ie proc6d6 
selon I'invention. et agir comma « pr^conditionnement » & la correction r6alls6e 
par la d6temilnation du tauxde couverture d'illumination pr6c6demment d6crite 
20 La fonctlon d'incr^mentatlon peut alo.^ §tre pond6r§e par ce terme de 

pr^condltionnement. de s^rte que Ie taux /(ic,^,..F) de couverture d'illuminatiori 
relatif S une pairs (source s . recepteur r ) est incr6mente par : 
i(td, tri s,r).M^x;s,r). 

L-evaluation du taux de couverture d'illumination dans Ie domaine migre 
25 selon l invention est optimum puisque applicable a n'importe quelle geometrie 
d'acquisition. 

Base sur une approche Inverse, Ie taux de couverture d'illumination l(^ p) 
estime par I'invention Concorde avec les r6sultats theoriques etablis par Bleistein 
Cela n'est pas Ie cas par les autres methodes proposees afin de resoudre un 
30 problfeme similaire. 
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Comme cela a d6ja 6t6 expos6 prec6demment, revaluation du taux de 
couverture d'illuminatlon propose par I'invention utilise seulement les quantit^s 
basiques disponlbles pour la migration classlquement r§alisee par sommation de 
Kirchhoff. 

5 De manldre avantageuse, revaluation du taux de cx)uverture d'illuminatlon 

selon I'invention peut etre r§alis6e tout en effectuant la migration des traces 
sismiques. 

Cela ri'est pas le cas de certaines des autres m^thodes proposees 
(Schneider et al. ; Laurain et Vinje) pour lesquelles ladite evaluation est realisee 

10 post-migration. 

Selon un autre aspect, I'invention propose egalement un precede de 
correction des amplitudes des donn^es sismiques enregistr6es au cours d'une 
prospection sismique 3D afin de compenser de I'effet d'une illumination non. 
uniforme des r6flecteurs du sous-sol, caract§ris6 en ce qu'il comporte les etapes 

15 consistant ii : 

. estimer le taux de couverture d'illuminatlon l(3c,p ) par le proc§d6 selon le 
premier aspect de I'invention d6taill6 d-dessus, 

• utiliser I'inverse \'^{x,p) dudit taux comme facteur de pond6ratioh d 
appliquer S chacune desdites amplitudes des donnees sismiques. >' 

20 Selon encore un autre aspect, I'invention propose egalement un proced§ 

pour la selection, en fonction de la cible d'une prospection sismique 3D, d'urje 
geometrie d'acquisition P|fmi une plurality de geometries d'acquisition, caracterise 
en ce qu'il comporte les etapes consistant § : 

• determiner, pour cliacune des geometries d'acquisition conslderees, le taux 
25 de couverture d'illuminatlon \{x,p ) estime conformement au,premier aspect 

de I'invention detailie cl-dessus, 
. selectionner la geometrie d'acquisition procurant le taux de couverture 
d'illumination optimum en fonction de la clble 

Enfin, on notera que le taux de couverture d'illumination \{x,p) estime 
30 confomiement au premier aspect de I'invention detailie ci-dessus, peut aussi etre 
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'r'*'"®® 3P''^s (a migration, pour 6liminer les empreintes d'acquisitipn. ainsi.que pouc 
deconvoluer I'image migree. 
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REVENDiCATIONS 

. Proc§d6 d'estimation du taux \{x,p) de couverture d'lllumlnation sismique 
dans le domaine 3D migre en au moins un point Image x , pour au moins un 
pendage de vecteur p, 

caract6ris6 en ce que, pour chaque paire (source s , r6cepteur r) du relev6 
sismique, on evalue le taux I(x,p;s,f) de couverture d'illumination relatif S 

ladite paire (source 2 , r^cepteur r ) par mise en oeuvre des §tapes de : 

• d6tenmlnation du temps de trajet de reflexion i,{xr{p) ', s, r) depuis 
la source s jusqu'au point de reflexion speculaire xr sur le reflecteur 
plan passant par le point image x et perpendiculaire au vecteur p 
de pendage. puis de retour vers le r6flecteur r , 

k partir du temps de trajet de diffraction td(x ; s, r) depuis la source 
s jusqu'audit point image jc, puis de retour jusqu'au r6flecteur r ; 

• d'incr6mentation dudit taux I{x,p\s,f) de couverture d'illumination 
relatif a ladite paire (source J, r6cepteur r) en fonction de la 
difference entre le temps de trajet de diffraction td(x; s, et |e 
temps de trajet de reflexion tr( xr ( ) ; s, r). 

2. Proc6d6 selon la revendication pr§c6dente, caract^rise en ce qu'il 
comporte en outre une etape de sommation de ciiacuri desdits taux 
/(jc,p;l,F) de couverture d'illumination relatif ^ une paire (source "s, 
recepteur r), de manifere a determiner le taux de couverture d'illumination 
total /(x, p) = Y,^ix,p;sJ) . 



1er cl§p6t 



20 



10 



15 



20 



3. Proced6 selon Tune des revendications pre.c§dentes,..caracteris6 en ce 
que. au cours de r6tape d'incrementatlon. le taux I(x,p;s,r) de couverture 
d'illuminatlon est incr6ment6 en utilisant une fonctlon d'inci^mentation i(td. 
V; S,r), selon Kx,p)^\{x,p) + l(td. tr; J,F). ladlte fonction d'lncr6mentatlon 
prenant en compte la dlff§rence entre le temps de trajet de diffraction U{x; 
s, r)et\e temps de trajet de reflexion U.Xr{p),s,r). 

4. Proc6d§ selon la revendlcation pr§c§dente, caract6ris6 en ce que la 
fonction d'incr6mentation i est fonction de I'ondelette sismique s(t). 

5. Precede selon la revendlcation pr6c§dente. caract§ris6 en ce que la 
fonction d'incr^mentation i s'exprlme en fonction de I'ondelette sismique s(t) 
selon : 

i(td. tr; s,r) - s{t^(x;s,r)-t^(x^(py,s,ry) . 

6. Proced6 selon la revendication 4. caract^rise en ce que la fonction 
d'incr§mentation i s'exprime en fonction de la derives s(t) selon le temps 
de I'ondelette sismique s(t) selon : 

Ktd.tr; J,F) = s(t^ix;s,r)-t,{x^(jp);s,ry). 



7. Proc6d6 selon I'une des revendications 3^6. dans lequel la migration mef 
en jeu une correction a priori w(xis,n du taux de couverture d'illumination. 
caract§ris6 en ce que le taux I(x,p;s,F) de couverture d'illuminatlon relatif 

a une paire (source s, r6cepteur r) est incr6ment« de i(td, tr; 
25 s,P).w(x;s,r). 
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8. Proc6d6 selon l*une des revendications pr6cedentes, caract^rise en ce 
que Tetape de determination comprend le d6veloppement en serie de 

Taylor, du second ordre, du temps de trajet de diffraction td{x; r) autour 

5 du point image x : 

9. Procede selon la revendlcation pr§c6dente, caract6rise en ce que le point 
de reflexion speculaire x^(p) est determine le long dudit reflecteur de sorte 

10 que le temps de trajet de diffraction audit point reflexion spSculaire x^(Jp) 

soit stationnaire, selon : a (v^^^ {x;s, r) + A^ j^ (3c; ?, r).{x^ Cp) - ^)) = 0 . 

10. Procede selon les deux revendications precedentes, caracterise en ce que 
le point de reflexion sp§culaire Xr et le temps de trajet de reflexion 

15 tr(xr ( ) ; 5 , r ) sont 6valu6s selon les expressions suivantes: 

x,{p)= x-M^F'^b; 

dans lesquelles, 

• M, une matrice^S x 2) decrite par deux vecteurs s'6tendant le long du' 
20 refelcteur, et done perpendiculaires au vecteur p de pendage ; 

• ^ , un vecteur (2x1) des derives du premier ordre du temps de trajet 
de diffraction le long du plan de reflexion : b = M ^ -{^Jd ) 

• F est une matrice (2 x 2) des derives du second ordre du temps de 
trajet de diffraction je long du plan de reflexion : F = M^.(A,^^^).M. 
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11. Proc6de selon I'une des revendications pr§cedentes, caracterise en ce 
que l'§tape de determination utilise les cartes des isochrones de migration 
Uix; s, r) specifiees pour toute paire (source s, r§cepteur r) impliqu^e 

5 dans la migration en tout point image x du domaine 3D migr§. 

12. Proc6de selon I'une des revendications pr6cedentes, caracteris6 en ce 

que I'estimation du taux \{x,p ) de couverture d'illumination sismique dans 

le domaine 3D migre est realist lors de la migration par sommation de 

10 Kirchhoff des donn6es sismiques enregistr6es au cours de la prospection 

sismique 3D. 

13. Proc6d§ de con^ction des amplitudes des donn6es sismiques 
enregistr6es au cours d'une prospection sismique 3D afin de compenser de 

15 I'effet d'une illumination non unifomie des r§flecteurs du sous-sol, 

caract6ris6 en ce qu'il comporte les etapes consistant d : 

• estimer le taux de couverture d'illumination \{x,p) par le precede 
selon I'une des revendications 1^12; 

• utiliser ['inverse \'\x,p) dudit taux comme facteur de ponderation d 
20 appliquer ^ chacune desdites amplitudes des donn6es sismiques. 

14. Precede pour la "selection, en fonction de la clble d'une prospection 
sismique 3D, d'une g6om§trie d'acquisition parmi une pluralit6 de 
geometries d'acquisition. caracterise en ce qu'il comporte les etapes 

25 consistant : 

• detemiiner, pour chacune des geometries d'acquisition consider^es. 
le taux de couverture d'illumination \{x,p ) par le proced6 selon I'une 
des revendications 1 a 12 ; 

• seiectionner la g§ometrie d'acquisition procurant le taux de 
30 couverture d'illumination optimum en fonction de la cible. 
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Pour chacune des paires (source s, r§cepteur r) impliquees dans la migration 



Pour tout point image x dans le domalne migr^ 



Pour toute composante de pondage "p s6lectionn6e pour le 
calculde I'lllumlnation I(x,p;s,r) 



D§tenmination du temps de trajet de reflexion Uxr{p) : s, r) 
au point de r6fle)don sp§culaire xr sur le r6flecte0r plan _^ 
passant par le point image x et perpendiculaire au vecteur p 
de pendage 



Incrementation du taux I(x,p;s,r) de couverture d'illumination 

>elatif a la pairs (source s . r6cepteur r ) en fonction de la 

difference entre : 

le temps de trajet de diffraction td(x ; 5 , r) et 
le temps de trajet de reflexion tr(xr(p ): s.r). 
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Pour chacune des paires (source?, recepteur r") impKquees dans la migration 



Pour tout point image x dans fe domaine migr4 



Pour toute composante de pendagep setectionnee pour le 
calcut de ritlumination /(x, p, r ) 



Determination du temps de trajet de reflexion tr (xrip); r ) 
au point de reflexion speculaire xr sur le reflecteur plan 
passant par le point image 3c et perpendiculaire au vecteur p 
de pendage 



Incrementation du taux/Jx,p, T, r ) de couverture d'illumination 
relatif a la paire (source?, recepteur r') en fonction de la 
difference entre : 

le temps de trajet de diffraction td(x; 7,T)et 
le temps de trajet de reflexion trixr (p);Xr). 
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